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ABSTRACT 

Knowing the masses of the components of binary systems is very useful for constraining the possible scenarios that could lead to their 
existence. While it is sometimes possible to determine the mass of the primary star, it is challenging to obtain good mass estimates 
of the secondary of a single-line spectroscopic binary. If the sample of such binaries is large enough, however, it is possible to use 
statistical methods to determine the mass-ratio distribution, and thus, the mass distribution of the secondary. Recently, the mass 
distribution of companions to extremely low-mass white dwarfs was studied using a sample of binaries from the ELM WD Survey. 
I reanalyse the same sample with two different methods: in the first one, I assume some functional form for the mass distribution, 
while in the second, I apply an inversion method. I show that the resulting companion-mass distribution can be as well approximated 
by either a uniform or a Gaussian distribution. The mass-ratio distribution derived from the inversion method without assuming any 
a priori functional form shows some additional fine-grain structure, although, given the small sample, it is difficult to claim that this 
structure is statistically significant. I conclude that it is not yet possible to fully constrain the distribution of the mass of the companions 
to extremely low-mass white dwarfs, although it appears that the probability to have a neutron star in one of the systems is indeed 
very low. 
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1. Introduction 


Extremely low-mass white dwarfs (ELMs) are thought to be 
the end products of binary star interactions and generally have 
companion (e.g., [Marsh, Dhillon & Duck 1995[ [Rebassa- 


Mansergas et al.||20H| l. Although still far from understood, bi¬ 
nary models predict that the companions to these ELMs are 
white dwarfs because the system has undergone one or two 
common-envelope phases (|Nelemans et al.|200T van der Sluys| 
et al.||20()6{ Woods et al.||2012' Toonen et al. 2014| l. Knowing 
the distribution of the companion’s mass to these ELMs could 
provide useful constraints on the various parameters that enter 
binary star evolution and the common-envelope phase, as well 
as to predict whether and when they will merge and, more gen¬ 
erally, how they will evolve ( |Dan et al.|2012[pCilic et al.|2012| l. 
Thanks to the ELM Survey, which has identified 61 ELMs and 


provided orbital parameters for 54 of therrj^ (see Brown et al. 
2013| and references therein), it is now possible to begin envis 


aging conducting a statistical analysis as was done by [Andrei 


et al.|P0i4) l. The problem is that ELMs are single-lined spectro¬ 
scopic binaries and as such, it is not possible to obtain the mass 
ratio directly (see, e.g., Boffin|20T0 2012 Cure et al.|2014 l, but 
one needs to apply statistical methods to derive the mass-ratio (or 
companion mass) distribution. [Andrews et al.| ( |2014[ ) developed 
a Bayesian probabilistic model to infer the companion mass dis¬ 
tribution for the above-mentioned sample, assuming a functional 
form - a two-component Gaussian, with one component repre- 


Andrews et al. 


j2014| used the 55 systems of 


Gianninas et al. 


(20141, 


which are essentially based on the list of 54 systems from jBrown et a .[ 
( [2013^ . Here I used these 54 latter systems because they come from a 
homogeneous sample. 


senting white dwarfs with masses between 0.2 and 1.44 Mq, and 
the other neutron stars with masses centred around 1.4 Mq and 
a standard deviation of 0.05 Mq. Using a Markov chain Monte 
Carlo algorithm, they found that their best fit is given by a pop¬ 
ulation of white dwarfs centred around 0.74 Mq and a standard 
deviation of 0.24 Mq, without a neutron star. This is quite an in¬ 
teresting result that also indicates that in contrast to population 
synthesis models, the majority of companions to ELMs are CO- 
core WDs, and not another He WD. As such, it is important to 
examine whether this result holds when using different methods, 
including when the functional form is not fixed a priori. This is 
what I present here. 


2. Methods 

The derivation of the orbital elements for a single-lined spec¬ 
troscopic binary (period, radial velocity amplitude, and eccen¬ 
tricity) allows obtaining the spectroscopic mass function, f(m), 
which is a combination of the masses of the two components and 
the (unknown) inclination of the orbit on the line of sight, i: 


f(m) = 


Ml 

(Ml M2)2 


sin^ i. 


where Mi is the mass of the primary (in this case, the ELM), and 
M 2 , the mass of its companion. If Mi is known, as it is the case 
here, then one can rewrite this as a function of the mass ratio, 
q — M 2 /M 1 : 

Ml (!+#■ 
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The distribution of the logarithm of f{m) - or F - can be used 
to determine the distribution of M 2 or q. This was done here 
using two different methods, in which I always assumed that the 
inclination i is randomly distributed on the sky, that is, P{i) - 
sin i. I used the sample of Brown et al. ( 2013[ similar to that used 
by Andrews et al. ( |2014| l, but see the footnote), which provides 
a list of 54 systems with known f(m) and Mi. 

In the first method, I assumed a functional form for the distri¬ 
bution of M 2 : this is either a Gaussian with mean p and standard 
deviation cr, or a uniform distribution, defined between a lower, 
M 2 ,i, and an upper, M 2 ,u value of the companion mass. I then 
applied a Kolmogorov-Smirnov statistical test. I assumed a sin¬ 
gle population of companions, without distinguishing between 
a neutron star (NS) and a white dwarf companion population, 
for instance. This is based on the fact that |Andrews et al!]p014| l 
found the NS fraction to be very low, as I confirm here as well. 

In each case, I ran 10,000 Monte Carlo simulations, where 
Ml is distributed according to the observed distribution (see the 
online Fig. [^, M 2 according to the chosen distribution with one 
set of parameters, and i is assumed to be randomly distributed 
on the sky. For each sample of simulations, I calculated the cu¬ 
mulative distribution of log f{m), which I compared with the ob¬ 
served one. To do this, I calculated the largest deviation between 
the two distributions, D* - Z) 54 _ioooo- Running 10,000 simula¬ 
tions provides very good precision, and there would be no gain 
in running more, as the estimator D* saturates for large num- 
ber^ 

This estimator allows determining the probability with which 
the simulated and the observed distribution are extracted from 
the same population. Thus, if D* > 0.2628, the chance is only 
0.1% that the two populations are extracted from the same pop¬ 
ulation, and we may most likely ignore such solutions. For 
D* > 0.1274 and D* > 0.0506, these probabilities become 33% 
and 99.9%, respectively. The first value provides a 1-cr estimate 
of the parameters that are allowed, while the latter number can 
be used to estimate the 0.1% confidence interval of parameters 
that provide a very good match to the observed distribution. 

In the second method, I took a more direct approach and in¬ 
verted the distribution of Y to derive the mass-ratio distribution. 


3. Results 

3.1. Functional form 

3.1.1. Gaussian 

As mentioned above, I here assumed that the distribution of the 
companion mass, <[)(M 2 ), is given by a Gaussian: 

Mo — u 

< 1 )(M 2 ) oc exp(- 

and I determined for a plane 0 < p < 2,0 < cr < 1 the value 
of D*. The relation of D* with p is shown in the online Fig. 
where I also show the lines corresponding to D* = 0.0506, 
0.1274, and 0.2628. This figure shows that one can find simu¬ 
lations spanning the whole range 0 < p < 1.22 Mq that would 
lead to D* < 0.2628, while the l-cr range covers the values 
0.3 < p < 0.92 Mq . This indicates that it is difficult to constrain 
the parameters of such a functional form based on the small ob¬ 
served sample. It is clear, however, that the simulations with p 
between 0.7 and 0.8 Mq correspond to the lowest values of D*, 
with the minimum being around p-O.lb Mq and cr-0.21 Mq, 


^ This is because D„ „< oc oc ^ - for n' » n and n » 1. 

y nn' y n 
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Fig. 1. Standard deviation, cr, versus the mean mass of the companion, 
p, for all simulated samples where the companion mass is distributed 
according to a Gaussian, and which have D* < 0.1274 (in black) and 
those that have D* < 0.0506 (in green). The heavy, red dot shows the 
location p = 0.74, cr = 0.24. 



Minimum mass of companion 


Fig. 2. Highest mass of the companion as a function of the lowest mass 
of the companions for all simulated uniformly distributed samples that 
have D* < 0.1274 (in black) and those that have D* < 0.0506 (in green). 
The two masses are clearly correlated, and I show a linear fit to the green 
dots as the heavy, solid red line. 


leading to D* - 0.027. This value indicates that this simula¬ 
tion and the observed distribution are consistent at a very high 
level because the null hypothesis that the two populations are 
drawn from the same population can only be rejected at a level 
of 1.3 X 10-'2%. 

Of course, the standard deviation cr is correlated with p, 
and Fig. [T] shows the solutions that lead to D* < 0.0506 and 
D* < 0.1274, illustrating the allowed range. For D* < 0.1274, 
I found that the mean value of the companion mass is = 
0.70 + 0.12 Mq, while the mean value of the standard deviation 
is (Tm - 0.45 + 0.22 Mq. For D* < 0.0506, these values become 
Pm = 0.76 + 0.03 Mq and cr„ = 0.27 + 0.06 Mq. 

A analysis that computes the deviation of the computed 
and observed distribution of log f{m) provides a similar result, 
with the best fit being given hy p - 0.76 + 0.02 Mq and cr = 
0.27 + 0.02 Mq (see the online Fig.|^. 
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Fig. 3. Comparison between the observed distribution of the logarithm 
of the spectroscopic mass function (solid black line) and the best fits 
for the uniformly distributed (red dotted line connected by heavy dots) 
and Gaussian-distributed (green dashed line connected by open squares) 
companion masses. For the former, I used 3 / 2,1 = 0.25 Mq and M 2 ,u = 
1.28 Mq, while for the latter, I used p = 0.76 Mq and cr = 0.27 Mq. The 
top panel shows the fraction of systems, while the bottom panel is the 
cumulative fraction of systems. It is clear that both samples are good 
fits, given the intrinsic errors of the observed distribution. 


3.1.2. Uniform distribution 

I repeated the same analysis but using a uniform distribution of 
the companion mass, between a lowest (M 2 . 1 ) and a highest value 
(AT 2 ,u), which were assumed to be in the range 0 < M 2 ,i < 0.9 Mq 
and M 2 ,i < M 2 ,u < M 2 ,\ + 1.5 Mq. The results are shown in the 
online Fig. [7] and in Fig.[^ Again, the range of allowed values is 
very wide. If this is restricted to D* < 0.2628, the whole range of 
M 2.1 is allowed, while for M 2 ,u, it is restricted to values between 
0.5 and 1.95 Mq. The range becomes narrower for lower values 
of D*, as shown in Fig. 2, which also shows that the acceptable 
values of M 2 ) and M 2 ,u are correlated. A linear fit gives 

M 2 .U = -1.25M2,i + 1.60. 

The mean values of M 2 ) and M 2 ,u are 0.36 ±0.18 Mq and 
1.20 + 0.23 Mo for D* <’ 0.1274,’and 0.33 + 0.07 Mq and 
1.19 + 0.10 Mo for D* < 0.0506. The lowest value of D* 
(0.02793) is reached for M 2 ) - 0.25 Mo and M 2 ,u = 1.28 Mo. 
This value corresponds to a probability that the two distributions 
are not extracted from the same population of 4 x 10 '°%. Per¬ 
haps most importantly, this shows that a uniform distribution of 
the companion mass fits the data very well and cannot be dis¬ 
carded. With such a functional form, a Gaussian or a uniform 
distribution cannot be preferred; they do, of course, quite over¬ 
lap (see the online Fig. |^. 

This is further illustrated in Fig.[^ which shows the distribu¬ 
tion of log/(m) for the observed sample as well as for the two 
best fits of the functional form, a Gaussian, and a uniform distri¬ 
bution. It is clear that both simulated distributions are good fits 


to the observed distribution, while it is very hard to distinguish 
the results from the two different functional forms. 


3.1.3. Inversion method 


Using a functional form for the mass distribution allows a better 
control of the systematics of a method, but at the cost of risk¬ 
ing missing some interesting deviations. This is for example the 
case in the paper by |Andrews et al.| ( |2014[ ) in their test 4, where 
they compare the distribution they obtain for a sample of post¬ 
common envelope systems to that obtained from spectroscopy. 
Although they clearly reproduce the bulk distribution, they miss 
the tail and other details of the distribution (see their Fig. 3). 
Moreover, [Andrews et al. ( |2014| l also pointed out that for the 
sample of ELMs, their result “could indicate that the true WD 
distribution may not be exactly Gaussian” - and indeed the pre¬ 
vious section has shown that a uniform distribution also provides 
a good fit to the observations. It is therefore useful to consider 
exploring methods that do not require the a priori input of a func¬ 
tional form. This is the case of the Richardson-Lucy (R-L) inver¬ 
sion method, as used by [Boffin, Cerf, & Paulus|(|1993]l; [Boffin] 
( |20T0l[20T2t , and [Cure et al.[ ( [2014| l. 



Fig. 4. Distribution of the mass ratios as determined by the Richardson- 
Lucy algorithm (solid black line) and as obtained for the best fits with a 
Gaussian (red dotted line) and a uniform (blue dashed line) distribution 
of companions. 


I refer to these papers and references therein for a full dis¬ 
cussion of the method, and in particular to Cerf & Boffin ( 1994| l 
for a more formal presentation. Here, I just mention that the 
Richardson-Lucy method relies on the Bayes theorem on con¬ 
ditional probabilities and solves the Fredholm integral equation 
that links Y with q by an iterative scheme. It is important to note 
that we do not directly have the distribution of the companion 
mass, but only of the mass ratio. This is, however, a very impor¬ 
tant parameter for binary evolution models, and given that the 
mass of the primary is very peaked at 0.17 Mq, the distribution of 
the mass ratio can give an idea of the distribution of the compan¬ 
ion mass. The outcome of this method is shown in Fig.|^ where 
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I also compare it with the distribution derived from the best fits 
of the functional form method. For the latter, 1 derived the mass- 
ratio distribution by using the functional form for the companion 
mass and the primary mass distribution as determined by|Brown| 
[erar] ( |20T3] l. 

All three methods provide a rather similar mass-ratio distri¬ 
bution, with some small differences. The outcome of the R-L 
method gives a broad Gaussian distribution, but with two more 
pronounced peaks, around q ~ 2-2.5 and around q ~ 4-5. 
This may indicate that a simple functional form may be missing 
on some small structure in the data, although a much larger sam¬ 
ple is needed to be able to confirm these ( Boffin|[2012[ ). A K-S 
test indicates that the hypothesis can be rejected that the func¬ 
tional forms and the outcome of the R-L method are drawn from 
the same population at the level of 92% - a rather high number, 
but perhaps not convincing enough. Indeed, given the small sam¬ 
ple, this is at most a two-sigma result. However, it shows that the 
true mass ratio distribution (and companion mass distribution) 
may have a more complicated structure than any simple func¬ 
tional form we can think of. Only with much larger samples will 
be able to know this. 

At the suggestion of the referee, I have also examined the 
mass-ratio distribution derived with the R-L method when limit¬ 
ing to the least massive primary stars, that is, the 33 systems with 
masses lower than 0.2 Mq. The resulting distribution is shown in 
the online Fig. It shows a single-peaked distribution centred 
around q ~ A - 5, that is (given that we now examine systems 
with a primary mass Mi = 0.17 Mg), M 2 = 0.68 - 0.85 Mg. 
In the figure, the corresponding mass-ratio distributions for the 
functional forms were computed with a single value of the pri¬ 
mary mass. The outcome of the R-L method apparently agrees 
better with the Gaussian distribution, but given that we have now 
an even smaller sample, the data should not be overinterpreted 
because the three distributions are again compatible within 2-cr. 
This peak may correspond to the similar peak seen in the mass- 
ratio distribution seen for the entire sample, while the fact that 
there is a possible second peak in Fig. at a lower mass ratio is 
most probably linked to the more massive primaries. 

The very small excess of high mass ratios seen in Fig. |^and 
in the online Fig. [^should not be given too much importance: 
hrstly, it is a well-known effect of the R-L method to smooth 
the distribution, the effect becoming weaker depending on the 
number of iterations ( Cerf & Boffin|p^94) l, and secondly, their 
value is compatible with zero, given the size of the sample. 


4. Discussion and conciusions 


The results of Andrews et al. ( |2014 1 in the study of ELM WDs 
are very important, therefore I have reanalysed the same sam¬ 
ple of spectroscopic binaries they used with two different tech¬ 
niques. In the first one, I assumed functional forms and applied a 
K-S statistics. The obtained results were conhrmed by a;^f^ sta¬ 
tistical test. The parameters for the Gaussian functional form are 
very similar to those found using a different method by |Andre^ 
|et al.| ( |2OT4| ), but I showed that a uniform distribution of the com¬ 
panion mass can provide as good a fit to the observed data and 
that the range of parameters allowed is rather large, making it 
hard to provide a definitive answer as to the real distribution. If 
the uniform distribution illustrated in Fig. is more representa¬ 
tive of the companion mass distribution, then more double He 
WD binary systems may be expected, such as the eclipsing dou¬ 
ble white dwarf binary CSS 41177 ( |Bours et al.|2014| l: the uni¬ 
form distribution shown in this figure leads to a 24% probability 
(this is the fraction of systems that have a secondary mass lower 
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than 0.5 Mg), compared to 16% for the Gaussian ht shown in 
the same figure or as derived by [Andrews et al. ( 2014| l. If I take 
into account all the possible values at the 1-cr level (i.e. those 
with D* < 0.1274), I also derive a 26% probability for both my 
Gaussian and uniform distributions. 

In addition, I applied an inversion method to derive the mass- 
ratio distribution, without the need to assume any functional 
form. The results are compatible with those derived from the 
functional form, although they seem to indicate some additional 
fine-grained structures. 

An important question that Andrews et al. ( |2014| l addressed 
is the possibility to have a neutron star as a companion to the 
ELM WD. [Andrews et al.|p014| l found that the probability of 
this is very low. This is also confirmed by my results. Eor exam¬ 
ple, the uniform distribution of companion masses indicates that 
the range allowed for D* < 0.0506 ends at 1.4 Mg. If this is re¬ 
laxed to D* < 0.1274, however, the highest companion mass can 
be up to 1.7 Mg, but all in all, the probability to have one system 
with such a high mass (> 1.44 Mg) is very low: for the uniform 
distribution, the 1-cr probability is zero, while for the Gaussian 
distribution it is 6.5%. 


Acknowledgements. I would like to thank the referee, JJ. Andrews, for a careful 
reading of the manuscript and for providing suggestions to improve the paper. 
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Fig. 5 . Distribution of primary mass in the sample of [Brown et al. 
(( 20 T 3 ). The sample contains many objects with primary mass equal to 
0.17 Mq. 
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Fig. 7. (top) Maximum deviation as a function of the highest mass of 
the companion if the latter is uniformly distributed between a lowest and 
highest mass, (bottom) The same, but as a function of the lowest mass. 
The lines are drawn at the same value of D" as in the online Fig.|^ 


Fig. 6. Maximum deviation between the observed sample with N = 
54 and a synthetic sample containing N' = 10,000 objects, with the 
companion mass being distributed as a Gaussian with a mean mass as 
given by the abscissa. The blue dotted line is drawn at D* = 0.2628: all 
simulations with a higher value of the estimator can be rejected at the 
99.9% level as being drawn from the same population as the observed 
distribution. The magenta dashed line is drawn at D* = 0.1274 and 
provide a l-tr estimate of the allowed parameter range, while the red 
solid line is drawn at D* = 0.0506. 
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Mean mass of companion 



Fig. 8. Reduced map in the mean companion mass vs. standard devi¬ 
ation for an assumed Gaussian functional form of the companion mass 
distribution - contours are shown at lev els of 1. 1.5, 2, and 2.5 . The red 
dot shows the position of the best fit of |Andrews et al. (2014 l The 
map confirms the result that was obtained with the K-S statistics. 


Fig. 10. Distribution of the mass ratios as determined by the 
Richardson-Lucy algorithm (solid black line) for all systems with a pri¬ 
mary mass below 0.2 Mq, and as obtained for the best fits with a Gaus¬ 
sian (red dotted line) and a uniform (blue dashed line) distribution of 
companions. 



Fig. 9. Comparison between the two best functional forms that lead 
to the lowest values of D*: a Gaussian distribution with p = 0.76 Mq 
and a = 0.27 Mq (solid black line) and a uniform distribution with 
M 2 ,i = 0.25 Mq and M 2 ,u = 1.28 Mq (dotted red line). 
























































